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Abstract: In this paper we discuss how the peculiar properties of twisted mass lattice 
QCD at maximal twist can be employed to set up a consistent computational scheme in 
which, despite the explicit breaking of chiral symmetry induced by the presence of the 
Wilson and mass terms in the action, it is possible to completely bypass the problem 
of wrong chirality and parity mixings in the computation of the CP-conserving matrix 
elements of the AS = 1,2 effective weak hamiltonian, and at the same time have a positive 
determinant for pairs of non-degenerate quarks, as well as absence of 0(a) discretization 
effects in on-shell quantities with no need of improving lattice action and operators. 
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1. Introduction 

In this paper we want to give an explicit set of rules that, within the computational scheme 
offered by the maximally twisted lattice regularization of QCD Q |j, allows to extract from 
simulation data the matrix elements of the CP-conserving part of the AS = 1,2 effective 
weak hamiltonian (7i^ =1 and 7i^ =2 ) between pseudo-scalar mesons with no contamina- 
tion from mixing with operators of wrong chirality or parity, nor 0(a) discretization errors. 
In contrast, a rather complicated mixing pattern has to be considered in order to construct 
multiplicatively renormalizable chirally covariant operators if, instead, the standard Wilson 
fermion regularization is employed Q. 

To achieve these remarkable results we will have to treat on different footing internal 
(sea) and boundary (valence) quark lines. Non-degenerate sea quarks are introduced in 
pairs and regularized as described in Q . In this way one gets a real, positive determinant 
with no vanishing eigenvalues for non-zero values of the quark mass, thus also solving 
the well known problems with the spectrum of the standard Dirac- Wilson operator. We 
immediately note that, through quark loops, all lattice quantities at finite o will obviously 
depend on the set of Wilson parameters, {r p }, associated to all the dynamical pairs of 
quarks (p = £(ight), /i(eavy), . . .) present in the action. 

As for valence quarks, a precise choice of their flavour structure and regularization 
scheme has to be made. In this paper we will discuss in detail the case in which valence 
quarks are introduced a la Osterwalder and Seiler Q , but one could equally well imagine 
to employ fermions obeying the Ginsparg- Wilson (GW) relation || for this purpose @. 
In this case 0(a) improvement and no wrong chirality mixing is also achieved, though 
at the expenses of a fairly larger computational burden compared to what is necessary 
for the inversion of the standard or twisted Wilson-Dirac fermion operator. It should be 
immediately observed, however, that the effort to invert, say, the overlap ||, ^] operator is 
not expected to be larger than that of stochastically evaluating the twisted determinant. 

In this paper we want to stick to Wilson-like fermions. We will show that one can 
still attain the ambitious goals mentioned above within the twisted mass regularization 
of QCD, with a judicious choice of the structure of the regularized valence quark action. 
More precisely, we propose to describe each valence flavour, /, with an action of the type 

S%£%f,qf,W = a^q f (x)h-V + e^W a (r f )+mf\ a f (x) 

d —~" L J LU=n/2 

x 

= fl4 Vf( x ) 7 • V - i^W„(r f ) + rrif q f (x) , (1.1) 



where 



7-V=iX)^(V^ + V^), (1.2) 



2 

Wcr(r/) = -a^E V ^ + M cr(r/), (1.3) 

with the Wilson parameter for the flavour / and M cr (rj) the corresponding critical 
quark mass. As remarked above, in the full theory M cr (rj) also depends on the set of sea 
quark parameters {r p }, though for short we have not indicated this dependence. 1 We will 
specify below more carefully what is to be meant by "critical mass" in the present context. 
We wrote the valence quark action in the so-called "physical basis" where my is taken 
to be real (and positive). 

The key feature we will exploit in the following is the freedom to regularize each 
flavour with an appropriately chosen value (of the sign) of the r parameter in front of the 
corresponding Wilson term. 

One might wonder at this point whether the scheme we are proposing is a consistent 
one, precisely in view of the fact that we are regularizing quarks in different manners 
according to whether they are sea or valence quarks. It is not difficult to convince oneself 
that, provided renormalized sea and valence quark masses of the same flavour are made to 
coincide, this is indeed so. The idea of the argument is to introduce compensating lattice 
ghost fields (a la Morel |jlC| 1) to exactly cancel the fermion determinant of the valence 
quarks and exploit the symmetries of the extended (local) theory |11] to prove that the 



Green functions of operators constructed in terms of only ordinary fields have the expected 
renormalization properties as the lattice spacing, a, is sent to zero (see appendix [A] for 
some detail). 

We will show that, by replicating some of the valence flavours (and accompanying 
ghosts) and regularizing them with Wilson terms of suitably chosen signs, the resulting 
lattice action turns out to have a sufficiently large symmetry so as to imply the cancellation 
of all the unwanted mixings in the lattice computation of the matrix elements of the CP- 
conserving AS = 1 and AS = 2 effective weak hamiltonian. It is important to remark 
that all these modifications of the lattice action can be made without spoiling the "magic" 
improvement properties of maximally twisted LQCD 0]. 

A preliminary study of the scaling properties of maximally twisted LQCD in the 



quenched approximation has been carried out in ref. [12], where measurements of the 
pion decay constant and the mass of the vector particle have been performed in the inter- 
val (5 = 5.85 ~r~ 6.2. Data extrapolate very smoothly to the continuum limit with residual 
discretization errors consistent with the expected 0(a 2 ) behaviour. 

Concluding, we wish to cite the nice paper of ref. Jl3| in which an interesting step in 



the direction of alleviating the renormalization problems of the Ti. eS ' operators within 
the framework of twisted mass LQCD (tm-LQCD) has been taken. The authors show that, 



1 Unless differently stated, we employ here the notations of ref. [0|. In particular, the ubiquitous depen- 
dence on the bare gauge coupling, g%, will always be omitted. 



with an appropriate choice of the twisting angles of the various flavours, in the computation 
of K — ► 7r transitions it is possible to remain with at most a linearly divergent counterterm, 
which in some partially quenched approximation can be even reduced to a finite one. 

The plan of the paper is as follows. In section || we discuss how valence quarks must 
be introduced and regularized in order not to loose 0(a) improvement and at the same 
time be ready to get cancellation of wrong chirality and parity mixings in the relevant 
physical amplitudes. In section |3] we present the strategy for the calculation of Bk, showing 
how in practice one can get rid of 0(a) discretization errors and completely dispose of 
the contributions potentially coming from the mixing of the bare operator H^ =2 with 
operators of unwanted chirality and parity. In sections || and || we extend this procedure 
to the calculation of the K — > tt and K — > tttt matrix elements of the CP-conserving 
^eff S=1 °P era t or - A few concluding remarks can be found in section |(| Technical details 
are relegated in some appendices. 



2. Sea and valence quarks 

In this paper we will consider QCD with four active quark flavours, a situation which rep- 
resents a rather realistic framework for the computation of the hadronic matrix elements of 
the effective weak hamiltonian. For the reasons explained in the Introduction we regularize 
in a different way sea and valence quarks and thus we take the lattice action to be of the 
form 



S = S 9 [U}+ E ^ /2) [^^^]+E[ 5 os /2) [^^^]+^ /2) ^/^] 



(2.1) 



P =i,h f 



where U are the gauge link variables and S g [U] stands for any sensible discretization of the 
pure gauge action. We briefly comment on the three "fermionic" bits of this action. 

• Sea quarks (p = £, h) are arranged in non-degenerate pairs: ipi with components 
u sea ,,dsea (^ight doublet) and iph with components s se a,c sea (heavy doublet). As pro- 
posed in ref. Q, they will be described by the maximally twisted action 2 

S tm [VV> $v> V) = fl4 E 7 • V - i^TiW CT {r p ) + m p - e p r 3 ip p (x) (2.2) 

X 

with n and T3 acting in flavour space. Integration over the sea quark fields leaves 
behind a real positive (for \e p \ < \m p \, p = £,h) determinant (see ref. Q and sec- 
tion U). 

• Valence quarks are denoted by qf. Each valence quark will be described by the OS 
lattice action, Sq£ \[qf,qf,U], specified in eq. (1.1). Since it will turn out to be 
necessary to introduce "replicas" of the usual u,d,s,c quark fields, we immediately 
notice that the range of variability of the index / will cover the more extended set of 
values u, u', it", . . . , d, d', d", . . . , s, s', s", . . . , c, c', c", .... 



2 With respect to ref. m we have inverted the sign in front of the mass splitting term. 



In order to have a local representation of correlators with insertion of valence quarks, 
each qf will be introduced with an accompanying ghost, cpf, which has the only 



purpose of exactly canceling the corresponding valence determinant [10|. The action 
of each ghost reads 

s£ /2) [</>/, U] = a 4 0/0*0 sign(m/) [7 • V - ij 5 W CI (r f ) + m f ] </, f (x) . (2.3) 

X 

The (complex) field (j)f is an euclidean commuting spinor in the fundamental repre- 
sentation of the gauge group with precisely the same value of the mass parameter, mj, 
and the same lattice regularization as the quark field, qf, it corresponds to. Since we 
will never need to consider correlators with insertions of ghost fields, nor graded field 
transformations (which turn valence ghosts into valence quarks and viceversa) , we do 
not need to introduce the sophisticated and powerful formalism developed in . 

The factor "sign (m/)" in the ghost action (|2.3| ) is inserted in order to guarantee the 
convergence of the gaussian integral over the 4>f field in case mj < [pd] ]. Indeed, with 
the purpose of discussing renormalizability and 0(a) improvement of the correlation 
functions computed with the action ( |2.1| ), in the following we will introduce spurionic 
symmetries that involve the inversion of the sign of mt. Thus the case of negative 



m 



f arises, even if one starts with mj > 



3 



As it is formulated, the lattice model with action ( |2.1| ) can be viewed as an euclidean 
statistical theory where the approach of correlation functions to the continuum limit can be 
discussed on the basis of standard power counting and symmetry arguments. In particular 
it turns out that renormalized correlation functions admit a well defined continuum limit 
provided, besides suitably rescaling the bare gauge coupling and setting the critical mass 
parameters of sea and valence quarks (and corresponding ghosts) to the appropriate values, 
for each flavour renormalized masses of sea and valence quarks (and ghosts) are given the 
same numerical value. Actually a continuum limit exists even if renormalized sea and 
valence quark masses of the same flavour are given different values: however, we will not 
consider the possibility of partial quenching in this paper. 

Energy levels and operator matrix elements of QCD with four quark flavours can be 
extracted from correlation functions of the model ( |2.1[ ) with neither ghost nor sea quark 
fields. As for correlators with sea quark fields, they are also perfectly well defined and are 
needed, for instance, to compute the renormalization constants Zp and Z$ (see eq. (^1]) 
below and appendix [A]). 

2.1 Symmetries and renormalization 

Since the lattice regularizations of sea and valence quarks are not identical, implementing 
the equality of sea and valence renormalized quark masses requires the knowledge of the 
relation between bare and renormalized masses. Such a relation follows, as usual, from the 



3 Notice that the factor sign(m/) does not introduce any spurious non-analytic behaviour of the ghost 
determinant in m/, as its presence only amounts to rescaling the determinant by an overall field-independent 
factor, which in any case drops out from (connected) correlators. 



analysis of the symmetries of the model described by the action (|2.1| ). The result of this 
analysis, which is presented in appendix |A[ can be summarized as follows. 



A unique (dimensionless) lattice function, / C r(n; r 2, r 3), exists in terms of which all 
critical masses can be expressed. Besides the coupling constant (which, as we said, 
we do not show explicitly), / cr depends on the Wilson parameter, r\, of the quark 
whose critical mass is being considered and on the Wilson parameters of sea quarks. 
The former — which appears only in the "external" lines of the diagrams relevant for 
the study of quark mass renormalization — can be rj, or either rg or r^, depending 
on whether we are dealing with valence or sea quark masses. The further dependence 
upon the Wilson parameters r2 = and = is the one that comes from "internal" 
loops. 

In the following for short the dependence on sea quark r-coefficients (r£ and rh) will 
always be understood. Thus the relevant formulae for the critical masses will be 
written in the form 

sea aM cr (r p ) = f cr (r p ; r E , r h ) , p = £,h, (2.4) 

valence aM cr (r f ) = f cr (r f ;r e ,r h ) , all f's . (2.5) 

It may be worth noting that also for valence quarks the value of the critical mass is 
determined by the requirement that qfy^W^qf (eq. (|1.3|) ) is a truly dimension five 
(irrelevant) operator. 

The renormalized quark mass parameters are related to their bare counterparts by 
the relations 



m p 



p = Z p l (r p )m p ±Z s 1 (r p )e v , p = £,h, (2.6) 
rhf = Z m (r f )m f , all f's , (2.7) 



where ± refers to the high and low mass components of each sea quark pair. Also in 
the above renormalization constants the dependence on the Wilson parameters that 
enter only through internal loops will be understood in the following. In other words, 
we will write Zj(r) as a shorthand for Zj(r; r£, r^), I = m, P, S, . . .. 

Zp(r p ) and Zs{r p ) are the renormalization constants of the (non-singlet) sea quark 
bilinears V> p 75T2 ) 3?/'p and i) p T2^ p -, respectively. They can be determined by requiring 



that the chiral Ward-Takahashi identities (WTI's) of the model (2.1) with insertions 
of sea quark field operators have the form expected in the corresponding continuum 
target theory. 

Z m {ff) is the inverse renormalization constant of the operator q/qj- The renormal- 
ized expression of the latter, in the model we are considering, is Z~ 1 (rj)[g/gj + 
mfCq q (r f ,ora/)a" 2 ], where the dimensionless coefficient Cq q is even in rj and am/. 4 



These properties of Cq q follow from the spurionic invariances (eq. (2.16)) and Vs x (M — > —M) 



(|A.9D and (2.12)) 



Z m (rf) is even in rf, as it follows, for instance, from the argument given below 
eq. (|2.21| ), and obeys the relation 

Z m (r) = Z m (r;r e ,r h ) = Z p l (r-r t ,r h ) . (2.8) 

Conservation of each of the N valence flavour quantum numbers is ensured by the 
invariance of Q2.1|) under the U(l) transformations (/ = !,..., N) 



x | q f (x) -» e ie fq f {x) , q f (x) ?/ (x)e ie f ^ g ^ 

with #y a real parameter. Larger vector valence flavour symmetries of course exist if 
several valence quark species have identical Wilson parameters (and bare masses). 

The flavour singlet axial current Js^ = J^/Li <Z/7a»759/; made out of exclusively va- 
lence quarks, exhibits the correct anomaly, proportional to JV Jl 



Generalizing the arguments given in ref. one can prove that / cr (^i; J"2 5 r 3) is odd in r±. 
For instance, if we take r\ = rf, proving that f cr is odd in its first argument is equivalent 
to prove that M cr (r/) is odd under rj — » — r/. This result follows (see appendix [A| for 
details) from the observation that under the changes of integration variables 



K 



5/ 



IX 



q f (x) -» ^(x) = 75^/(5 

-> q' f (x) = -q f (x) l5 (2.10) 
k <t>f(x) -> <fi' f ( x ) = 750/0*0 

the valence quark (plus ghost) action is invariant if we also invert the sign of rf and rrif (2J. 

An analogous argument applies if we were to consider the critical mass of sea quarks. 
In this case one has to perform the change of variables 

K 5 ■ { ^- p ^ ~* = l5 ^ p ^ (2 11) 

and exploit the fact that for each value of p the fermionic determinant is an even function 
of the corresponding Wilson parameter, r p (see section 2.2). 

In the following for any r G {re, r^; . . . , Tt , . . .} it is understood that W cr (r) is of the 
form (|l.3|) with M cr (r) the appropriate odd function of r specified in eq. fl2.4|) or eq. ([T^). 
We will thus label expectation values by the relevant Wilson parameters {re, r^; . . . , r/, . . .} 
and quark masses {me, eg, rrih, e^; . . . , mj, . . .}, for which we introduce the compact notation 

M = {rrif , all f's; m p , e p ,p = £, h} , R = {rf , all f's; r p ,p = I, h} . (2-12) 

2.2 Sea quark determinant 

Integration over sea quark fields gives rise to the product of determinants 

D = H B p , (2.13) 

p=e,h 

D 



v 



Det [7 • V - i75TiW cr (r p ) + m p - e p r 3 . (2.14) 



As proved in Q , each factor in D is even in m p and e p . From this result one can show that D p 
is also even in r p . To prove this property we follow again the strategy developed in ref. 0|. 



By performing in the functional integral that defines D p the change of variables (2.11), we 
get the identity 



Det 



7 • V - i-y 5 TiW cr (r p ) + m p - e p r 3 



Det 



7-V + ^75TiWcr(rp)-mp + e p T3 . (2.15) 



Since D p is even in m p and e p , while W CT (r p ) is odd in r p , we get the thesis. It then follows 
that all the correlators with insertion of operators made out of only valence quark and 
gluon fields, and thus all the quantities derived from them, are even in r p , m p and e p , 
p = £,h. In particular such parity properties imply that f C r( r i] r 2, r 3) is an even function 
of its last two arguments (r2 3 = r p ,p = £, h). 

2.3 Renormalization constants and O(a) improvement 



The symmetries of the action ( |2.1| ) relevant for analyzing the -R-dependence of renormal- 
ization constants and leading cut-off effects are 

• the spurionic IZ^ transformation acting on each flavour / separately 

Tlf f = K 5f x (r, -r f ) x (m f -m f ) , (2.16) 
with 1Z§f given by eq. ( 2.10| ); 

• the spurionic 7^5? transformation acting on the p = I and p = h sea quarks separately, 

= ^5 P x (r p -> -r p ) x (m p -> -m p ) x (e p -> -e p ) , (2.17) 
with lZs p given by eq. ( 2.11 ); 

• the overall transformation acting on all types of fields 



v d x Y[n 5f x [J ^5 P 

/ P=l,h 



(2.18) 



with 



V d : < 



ip p (x) 
q f (x) 



Ul(-x - afi) 
e m / 2 ip p (-x) . 
e 3i7T / 2 q f (-x) , 
e 3t7T / 2 cj) f (-x) . 



if> p (x) 
q f (x) 



/2 A 



= 3ivr/2 



q f (-x) 



p = £,h 
all fs 
all fs . 



(2.19) 



We also note the existence of the spurionic symmetry obtained as the product of all imparl- 
ances of the type ( 2.16 ) and ( 2.17 ), namely 

K L x (M —M) = Y[H 5f x Yl K 5p x(R^ -R) x (M -» — M) . (2.20) 



o 

to 
o 
o 

o 
-J 
o 



-8- 



Since the formal continuum action we are aiming at, of which ( |2.lD is the regularized version, 
will have to be invariant under the (discrete, non-anomalous) spurionic transformation 



K' 



cont 



x (M -» —M) 



l[n 5f x H K 5p x (M 

/ P=l,h 



-M) 



(2.21) 



it follows that in any renormalization scheme, besides Z g , the renormalization constants 
Zp, Zs and Z m , appearing in eqs. fl2.6|) and (|2.7l ), and those of all the local operators whose 
bare expression has definite transformation properties under 7£ L , will be even functions of 
R. This is, in fact, necessary for the renormalized operators to behave under the latter 
transformation as the corresponding bare ones. Furthermore, recalling that the fermionic 
determinant D p (eq. (2.14)) is even in r p , one concludes that all the above renormalization 
constants will be separately even in (or trivially independent of) r p , p = £,h and r/, all 
/'s. 

2.3.1 0(a) improvement via Wilson averaging 

By a direct generalization of the argument given in ref. the above symmetries and 
remarks imply that VKA's of expectation values of multiplicatively renormalizable (m.r.), 
multi-local, gauge invariant operators (with no ghost fields), 0(x±, X2, ■ ■ ■ , x n ), are 0(a) 
improved, i.e. 



(0(xi, . . .,x n )) 



WA 



(M) 



1 



(0(xi, . . .,x n )) 



(8(R)(0(x l ,...,x n )) 



(R,M) 
cont 



{-R,M) 



(M) 



+ 0(a z 



(2.22) 



where it is understood that the space-time arguments at non-zero relative 

distance in physical units. The multiplicative factor Co(^) ^ s a coefficient (even under 
R — > —R) whose presence is necessary to match the lattice to the corresponding continuum 
expectation value. 

Analogous results can be derived from eq. ( 2.22] ) for energies and operator matrix 
elements, namely 



£ M (k; R, M) + £ M (k; -R, M) = 2E%?(k, M) + 0(a 2 ) , 

(h,n,k\0\h',n',k')\ (RM) + 
+(a,n,k|0|/ l / ,n / ,k')| ( _ i j iM) = 2C0(-R)(/i,n,k|O|/t / ,n',k / ) + 0(a 



cont 



(M) 



(2.23) 



(2.24) 



where i?h )fl (k; R, M) is the energy of \h, n, k) | (R t M) j the n-th eigenstate of the lattice quan- 
tum hamiltonian (which we assume to exist) with three-momentum k and unbroken quan- 
tum numbers h, while O denotes an arbitrary m.r. and gauge invariant local operator. 

2.3.2 0(a) improved results from a single simulation 
The action ( p.l| ) admits the further spurionic symmetry 

V x (R -> -R) , (2.25) 



where V is the parity operation. The latter, with the definition xp = (— x, t), has in the 
present context the form 



V : < 



U (x) -» U (x P ) , U k (x) -» [/l(s P - afc) , k = 1,2,3 

tpp(x) 7o^ P (a;p) , ^ p (x) -> ip p (x P )j , p = £,h 

q f (x) 7og/(x P ) , qf{x P )^ , all /'s 

, </>/(x) -> 7o<^/(xp) , all /'s. 



(2.26) 



The spurionic symmetry ( |2.25| ) is particularly important, as its existence entails the pos- 
sibility of associating a definite parity label to m.r. local operators and to eigenstates of 
the lattice hamiltonian. 5 By arguments similar to those given in ref. [|| one can prove the 
validity of the symmetry relations 



E h>n {k; -R, M) = E h J-k; R, M) , 



(h,n,k\0\h',n',k h 



(—R,M) 



rjhnVor]h'n'{h,n, -k\0\h',n', -k' 



(R,M) 



(2.27) 
(2.28) 



where t)q, r]h n and denote the parity labels of the m.r. operator O and the external 
states, respectively. Inserting eqs. ( 2.27 ) and ( |2.28|) in ( [2.23|) and ( 2.24| ), one ends up with 
the formulae 



E hj7l (k; R, M) + £ M (-k; R, M) = 2E^(k, M) + O 
{h,n,k\0\h',n',k')\ {RM) + 



(2.29) 



+VhnWVh>n'{K n, -k\0\ti , n' , -k')\ {RM) = 2(g(R){h, n, k\0\h',n', k!) + 0(a 2 ), (2.30) 

v ' ' (M) 



cont 



which show that the WA's ( ^23|) and fl2.24j) can be evaluated in terms of correlators and 
derived quantities obtained from a single simulation with Wilson parameters R. 

2.4 Some conclusive observations 



It is worth observing that in the lattice model ( |2.1| ) the fermionic determinant ( 2.13| ) 
introduces neither lattice artifacts linear in a nor parity breaking effects. This can be seen 
by considering expectation values of pure gauge (multi-local, gauge invariant) operators, 
O pg . Such expectation values depend on the sea quark Wilson parameters, rg and r^, and 
are even functions of them, because D = T>£ is even, but are obviously independent of 
the valence quark Wilson parameters. We thus have 



\Opg\X\ , . . . , X n )) 



(R,M) 



{Opg (^1 , • • • , 3?n)) 



(-R,M) 



(2.31) 



from which, comparing with the W A formula ( 2.22 ) written for O — > O pg , we conclude that 
the expectation values of pure gauge operators trivially coincide with their Wilson average 



5 The assignment of a parity label to eigenstates of the hamiltonian can be done in close analogy with 
the procedure outlined in sections 5.2 and 5.3 of ref. p), where the case of twisted mass QCD with one pair 
of mass degenerate quarks was discussed. 



1— ' 

o 

O 
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and, hence, are automatically free of 0(a) cut-off effects. Moreover, if the operator O, 
has a definite parity, 770 = ±1, i.e. it satisfies the equation 



(O pg (xi,. . . 



x n )) 



(R,M) 



7]O pg (O pg (x 1P ,...,X n p)} 



(—R,M) 



(2.32) 



then for its expectation value, owing to eq. ( 2.31| ) 



one gets 



(O pg (xi, . . . 



x n )} 



(R,M) 



r)o pg {O pg (x 1P , ...,x nP )) 



(R,M) 



(2.33) 



which is precisely what one would expect in a parity-invariant theory. In other words, all 
parity breaking effects in correlators with insertions of valence quark operators (and not 
containing sea quarks or ghosts) come from the breaking of parity induced by the valence 
OS quarks themselves, and reduce to 0(a 2 ) in Wilson averaged quantities. 

The basic reason behind the validity of these properties is that, as a functional of the 
set of link variables, {U}, the (maximally twisted) fermion determinant, D = D [{£/}], is 
invariant under the replacement {U} — * {Up}, where {Up} is the parity transformed set 
of link variables. This result can be proved by performing in the functional integral that 
defines the determinant, D, the change of sea quark integration variables induced by the 
transformation in the second line of eq. ( |2.26| ) and recalling that D is even in r p ,p = £, h. 

All these observations may be of some importance in practice, because they suggest 
the interesting possibility of employing GW-type fermions to regularize valence quarks and 
maximally twisted Wilson fermions to describe sea quarks. 6 In this hybrid framework, 
whose computational cost does not look prohibitively large on the scale of today's avail- 
able computer resources, 0(a) improvement is automatic. Furthermore, the nice chiral 
properties of GW fermions are sufficient to ensure absence of all wrong chirality mixings 
in correlators with insertions of only valence quark and gluon fields. 

This last statement might not be obvious. Actually in the next three sections we will 
show that this remarkable result can be achieved even if valence quarks are regularized a 
la OS, provided some clever adjustment of the lattice regularization of (four-flavour) QCD 
and of its computational scheme is made. 

3. The Bk parameter 

The popular way to discuss the issue of CP violation in the Standard Model is through 
the study of the compatibility of experimental data with the constraints imposed by the 
unitarity of the CKM matrix (unitarity triangle) involving the third generation of quark 
flavours. 7 One of the key quantity entering this analysis is ex , which represents the amount 
of indirect CP-violation in K — > tttt decay. Exploiting its precise experimental value to 
constrain the non-trivial vertex of the unitarity triangle requires the knowledge of the 
parameter Bk- 

6 A very similar proposal, with ordinary (rather than maximally twisted) Wilson fermions playing the 
role of sea quarks, has been made in ref. [jjj. 

7 Good reviews on the subject can be found in ref. jl5| . For recent summaries of lattice results see the 
papers quoted in ref. [M. 



In the formal continuum euclidean QCD with four quark flavours (u, d, s and c), 
renormalized at the scale fj,, the kaon S-parameter, Bxif^), is defined by the formula 

(K°\6^fAA^)\K°) = ™MlFlB K {n) . (3.1) 
Bk can be extracted from the correlator (xq > 0, yo < 0) 

cSi\x,y) = ((d^s)(x)d^f AA (0)(d l5 s)(y)) , (3.2) 
where the bare expression of the renormalized operator Oyy= 2 AA is 

®VV+AA = {si^d){s^^d) + (s7 M 7 5 (i)(s7 M 7 5 cO (3.3) 

and the superscript (qcd4) denotes QCD with four quark species. 

We want to show that it is possible to set up a lattice regularization based on (twisted) 
Wilson fermions, where Bk can be computed in terms of a four-quark operator which does 
not mix with any other operator of the theory. 

In order to explain our strategy it is convenient to first consider at a formal level an 
auxiliary gauge model where, besides the four lightest sea quark flavours (u sea , d seSL , s sea 
and c sea ), each of the d and s valence quark species is duplicated. The valence quarks u, 
d, d', s, s' and c are introduced together with the corresponding ghosts, required to cancel 
the associated valence determinant. We will denote the correlators of this model by the 
superscript (4s6v ). It is not difficult to convince oneself that in this model the correlator 

C^QK^V) = ({^^{x)2Q$^ AA {0){d 7& 8){y)), (3.4) 

where 8 

Qvv+AA = {si^d){s^^d') + (s7 M 7 5 cO(s'7 M 75cO + 

+{siiid'){s'~i ll d) + (s7 m 7 5 gO(s'7m75cO , (3.5) 

contains the same physical information as the "target" correlator C^^ K \x,y) in four- 
flavour QCD. 

To see this we simply have to imagine to have regularized the two models in the same 
way, for instance, by using GW fermions for all quark fields. If the bare masses of the 
quarks of the same flavour are given the same value, i.e. if one sets 

(qcd4) (4s6v) (4s6v) (qcd4) Us6v) Us6v) Us6v) 

m y u = m y u = m^ sea , m\ = m y d = m K d , = , 

(qcd4) Us6v) Us6v) (qcd4) (4s6i>) (4s6i>) Us6v) fo „n 

ml = ml = mc sca , ml = ml = m\, = m y Ssc!l , (3.6) 

then by an elementary application of the Wick theorem one checks that the equality 

C { ^ Q i(x,y) = C^ d K \x,y) (3.7) 

holds (sea quark loops trivially contribute the same to the two sides of the relation) for 
any value of the lattice spacing. This result implies that the two operators, Oyy= 2 AA 

8 What we really mean with the superscript AS — 2 in Qvv+aa is As + As' = 2. 



in the (qcd4) model and Qyy^ 2 ^ in the (As6v) model, can be taken to have the same 
renormalization constants. The same is also true for the renormalization constants of the 
operator d^y^s in the (qcd4) model and the operators dj^s and d'^s' in the (4s6t> ) model. 
Moreover, it is obvious that the renormalized masses of the quarks of the same flavour are 
(can be chosen to be) equal in the two models. With this natural choice of renormalization 
conditions, the relation Q3.7| ) holds for the corresponding renormalized correlators and 
carries over to the continuum limit. 

The strategy we propose for the evaluation of Bk is based on eq. (3.7) and consists in 



considering the correlator C^q^{x } y) in the UV regularization provided by the maximally 



twisted action (|2.1|). The reason for selecting this UV regularization is that, with a careful 
choice of the Wilson parameters of valence quarks, Qyy^AA wm n °t with any other 
operator of the theory. 

3.1 A convenient lattice discretization of Ck'qk 

In view of the previous remarks, we are led to consider the lattice regularized gauge 
model ( |2,1D with valence quarks q u = u, q d = d, qd' = d', q s = s, q s > = s' and q c = c, 
each having an action of the form ( |l.ip (accompanied by the corresponding ghosts, (j) u , (fid, 
<f> d i, (p s , cf) s >, <j) c ) plus two sea quark doublets, ip £ = {u sea ,d sea ) and t/j h = (s sea ,c sea ). The 
details of the lattice formulation of this theory are as explained in general in section [| We 
will refer to this regularized theory as the {AsQv)\ m model. 

For reasons that will become soon clear we take the Wilson parameters of the valence 
down and strange fields to be related in the following way 9 

rd = r s = r d > = -TV . (3.8) 

Obviously, in order to obtain from the (4s6v )\ m model correlators that in the continuum 
limit are equal to those of four-flavour euclidean QCD, we must require the renormalized 
masses of the (sea and valence) quarks of the same flavour to take identical values, i.e. 10 

m~l = rh u , rhj = m d = rh d > , 

= m s = m s , , m\ = m c , (3.9) 

and be equal to those of the target theory, (qcd4). In terms of bare masses, while the 
conditions rhd = rhd' and m s = m s i are satisfied (owing to eqs. ( |2.7p and ( |3.8] ) and the fact 
that Z m (rf) is even in rj) by simply setting 

m d i = m d , m s , = m s , (3.10) 



the remaining four conditions in eq. ( |3.9[) imply a non-trivial one-to-one correspondence 
between m u , m s , m c and the sea quark bare mass parameters m?, €£, m^, e/j, as 
dictated by eqs. (O)-^). 



9 Other equivalent choices, like td — r s — —r d r = r s i are possible. All the arguments we give in this 
section can be immediately adapted to these other cases. 

10 Our notation is such that e.g. mf — m,d sca and — m SBCa . 



We now prove that in the lattice model {AsQv)\ m , with the set of quark masses, M, 
satisfying the conditions ( |3,9| ), the (bare) Bk parameter can be extracted from the corre- 
lator 

« 6 E4 4 S* m (^y) = « 6 E((^ S 0(^)2S^+5iA(O)(d75 S )(y)}| , (3.11) 

x,y x,y 



with the operator Qyv^AA given in eq. (3.5) and the (m.r.) pseudo-scalar (d'^s') and 
(J75S) densities 11 playing the role of interpolating fields for the K° and K° states. The 
proof is divided in two parts. We first show how Bk is related to the large time behaviour of 
the correlator ( p.ll ) and in the next section we prove that the operator (3.5) does not mix 



with any other operator with the same unbroken quantum numbers. A similar result was 



already derived within the standard tm-LQCD approach of ref. and exploited in ref. [17] 
in actual simulations. In the present scheme one gets two extra bonuses: an automatically 
0(a) improved determination of Bk (as it follows from the arguments we give at the end of 
this section) and a positive definite fermion determinant even in the mass non-degenerate 
case Q. 

On a lattice with very large temporal extension the leading contribution to the spectral 
representation of the correlator (3.11) in the {AsQv)\ m theory is given by 



-M K i\xo\-M K \yo\ 

AM k ,Mk ^l^5^1^ )(^ |2Q^ A |^ )(K |J75^|^) 



(3.12) 

(R,M) 



where K° (K'°) is the neutral kaon (anti-kaon) state created from the vacuum f2 by the 
operator —dj^s (s'j^d') and Mk {Mk 1 ) the corresponding meson mass. We remark that 
with the choice (|3.8[) of the (valence) Wilson parameters, the lattice kaon mass, Mk-, and 
kaon decay constant, Fk, defined in terms of the valence quark fields d and s may differ by 
0(a 2 ) effects from their counterparts, Mk 1 and Fk', analogously defined in terms of the 
twin valence quarks d! and s'. Keeping track of this difference, we write the relation which 
defines the bare Bk parameter in the form 

(K'°\ 2Q^f AA \K°) = ^-Mk'Fk'MkFkB^ . (3.13) 

The values of Mk, Fk, Mk 1 and Fk' can be obtained in the usual way from the study of 
the two-point correlation functions 

a 3 E((s7 5 o0(0)(d7075s)0z)>, a 3 E((^W)(0)(a ? 7o75^)(^)> • (3-14) 

X X 

The renormalized Bk parameter is finally given by 
with Zqas=2^ the renormalization constant of Qyy^ AA . 



BjcOi) = Z Q ^ AA { ai x)B h £™ , (3.15) 



11 



Mixing of (d'^/^s') with (d's 1 ) and of (d'yss) with (ds) is ruled out by the invariance of the lattice model 



under V5 x (M — » — M), where "Ps is defined in eq. (A. 9). More precisely these are O(a) effects that appear 



with coefficients proportional to a(m s + m c i) or a(m s — md) 



The correlation function ( 3,11 ) can be 0(a) improved via WA, according to the pre- 
scription ( 2,22| ), and the same holds for all masses and matrix elements that occur in its 
spectral representation. However, this is not necessary. In fact, thanks to the spurionic 
symmetry ( 2.25 ), Mk, Mk', Fk and F^ 1 (see eqs. ( [2.29 ) and ( 2.30| )), as well as all the 
quantities that can be extracted from the correlator fl3.11| ), are automatically free of 0(a) 
cut-off effects. This property follows from the fact that ( p. 11 ) is the expectation value of 
a m.r. operator with the insertion of two other operators taken at vanishing spatial three- 
momentum. Consequently, B^ re as defined above can be extracted from automatically 
0(a) improved data. The same is true for the renormalization constant of the operator 

)AS=2 
-VV+AA 



Qv^r i 2 /i a (see section 3.3 of ref. 



and section 2.3.2 of this paper). 



3.2 Proof of the absence of mixing in Qyy+'AA 

Absence of wrong chirality and parity mixing in the (renormalizable expression of the) op- 
erator Qvv+'aa ls P rove d relying on the symmetries of the lattice model ( |2.1[ ) supplemented 
by the specifications given in section 34. In particular the conditions ( |3,8D and ( 3.1 0| ) are 
assumed to be fulfilled. 

It is almost obvious that the operator Qyy^ 2 ^ can not mix with operators of dimension 
less than six, because of its flavour structure and the fact that the valence sector of the 



action (24) is diagonal in flavour. In particular, the valence flavour quantum numbers d, 
a", s and s' are all conserved, as it formally follows from the invariance of the action under 
the four U(l)'s vector transformations, If, f = d,d',s,s', of eq. fl2,9|) . We conclude that 
Qyv^AA can om y mrx with four fermion operators with As = As' = 1 and Ad = Ad' = — 1. 

In order to rule out mixing with other six-dimensional operators satisfying the above 
selection rules, one must use the following extra symmetries enjoyed by the model action 
we are considering, i.e. 

1. Ex(d, d') x {wid <-> rrid'), 

2. Ex5(s,s') x {m s «-> — m s i), 

3. C x [Ex(d, s) x (rrid «-> m s )] x [Exs(d',s') x (m<f «-> —m s >)], 

4. V 5 x (M -> —M), 

where C and V$ are defined in eqs. (A.l) and (A. 9) of appendix [A|, respectively, and 

Ex(/i,/ 2 ) : q h *-* q h q fl *-> q h (f) fl *-* (f) f2 , (3.16) 
' Qh -» 759/2 Qh -» -9/2 75 4>h -» 75<A/x 



Ex 5 (/i,/ 2 ) : 



Qf2 -» 759/! 9/2 



-9/i75 <t>h 75 ( 



The transformation Ex(/x,/2) simply represents the exchange of valence flavours f± and 
f'2 with its action appropriately extended to ghost fields. The transformation Exs(/i,/2) 
is nothing but the product 



Ex 5 (/i,/ 2 ) = Ex(/i,/ 2 ) x K 5fl x K 5h 



(3.18) 



Notice that, if rf x = Tf 2 , Ex(/i, /-j) is a symmetry of the action ( |2,1| ) when combined with 
(m/j — ► m/ 2 ). If instead = — r/ 2 , the action is invariant under Exs(/i,/2) x (m/j <-> 
-m/ 2 ). 

Symmetry 1 rules out the dimension six operators that, besides obeying the selection 
rules As = As' = —Ad = —Ad' = 1, are not even under ~Ex(d,d'). A basis of operators with 
the desired unbroken quantum numbers is provided by Qyyl^ AA itself and the operators [ 18 
19] that, with obvious notations, we list below 



OAS=2 f\AS=2 n A ^=2 /"}1CA 



f)AS=2 r> A S=2 r> A S=2 

^■VV-AA ' ^-SS±PP ' ^TT 

"TT 

For instance, we have (07^ = i [7^, 7^/2) 

Qva±av = {si^d){s'-i^d') ± (sj ll j 5 d)(s > j ll d') + 



+(s lfl d')(s' W d) ± {s W d'){s'^d) , (3.20) 

=2 

77 



Q I 2 = e^Ap (sa^ u d)(s'a X pd') + (sa^d 1 )^' a\ p d) . (3.21 



None of the (dimension six) operators in (3.19) can mix with Qyy+ AA thanks to the other 
three (spurionic) symmetries introduced above. In fact 

• symmetry 2 rules out QyyZ 2 AA and Q A g=p P (the latter are odd, while Qyy^AA ls 
even) ; 

• symmetry 3 rules out Q A j, =2 and Q A §^p P (the latter are odd, while Qyy^ AA ls 
even) ; 

• symmetry 4 rules out Qy A ^ 2 AV , Q A p± PS and Q A ~ =2 (the latter are odd, while 
Qvv+AA is even )- 

The symmetries we have employed to exclude mixing of Qyy^ AA with the operators ( 3.19] ) 
can be viewed as the analog of the usual (spurionic) CPS symmetry [jl9| translated in 
the present chirally twisted lattice framework. Indeed, if valence quark (and ghost) fields 
are rotated to the basis where the (subtracted) Wilson term appears in the standard (un- 
twisted) form, while the mass term is chirally twisted, the spurionic symmetry obtained by 
combining the transformations 3 and 4 above turns out to correspond to the CPS invariance 
of the standard Wilson theory. 

We conclude this section by observing that the whole argument about the absence of 
mixing in the operator entering the matrix element from which Bk is extracted could be 
repeated by duplicating only the valence strange quark, and not the down quark. If one 
does so, one also has to change the form of the operator in eq. (3.5) with the substitution 

Qvv+AA - i^d) (s' 7At d) + (s^j 5 d) (s 7 ^ 5 d) (3.22) 

and accordingly modify the expression of the interpolating field of K' from d'j^s' to d^s'. 
The reason why this is possible lies in the observation that, given eqs. (HD and fl3lcl ), 



the flavour d! is actually indistinguishable from d. 12 In section ^2 we have presented the 
line of arguments leading to absence of mixing in the apparently more complicated setting 
where also the down quark is replicated, because it turns out that in this way the proof 
goes through in a somewhat simpler fashion. 

4. The AI = 1/2 rule from K — > nn amplitudes 

A long standing puzzle in low energy particle physics, often referred to as "octet enhance- 
ment" or the AI = 1/2 rule ^0|, is represented by the surprisingly large experimental value 
of the ratio 

T(K — ► 7r7r)|Ar-i/o 

R(K - tttt) = -L- ~ 400 . (4.1) 

1 (K -> 7r7rJ| A/=3 / 2 

Although the rate of the AI = 3/2 K — > 7T7T weak decays can be reasonably well computed 
within our present understanding of field theory (based on O.P.E. and renormalization 
group arguments), theoretical estimates of the AI =1/2 amplitude seem to give a much 
too small value compared to the experimental number pl| . 

Due to the amount of seemingly well founded theoretical ideas involved in this calcula- 
tion, it is a crucial question to understand whether this discrepancy is due to our ignorance 
of non-perturbative QCD effects in the hadronic matrix elements of the composite quark 
operators appearing in the effective weak hamiltonian, or whether this problem should be 
rather taken as an indication of the existence of some kind of new physics waiting to be 
uncovered in kaon phenomenology. Lattice gauge theory represents the ideal framework in 
which this question should find an answer. There are, however, severe technical and con- 
ceptual difficulties in the process of establishing the proper strategy to accomplish this task. 

First of all, when working with euclidean metric, as one always does in lattice simula- 
tions, the procedure necessary to extract the kaon decay amplitudes of physical interest is 
complicated by IR subtleties arising from final (two-pion) state interactions [22]. To attack 



this problem new interesting ideas have been recently put forward [p3| , 24 1. Here we will 
not be concerned with such delicate issues. 

The main objective of this paper is rather to deal with the UV difficulties related to 
the construction of the renormalized effective weak hamiltonian operator on the lattice. 
We will show that tm-LQCD at maximal twist offers a computational scheme in which the 
problem of wrong chirality mixing, which afflicts the approach based on standard Wilson 
fermions |3], |25|], 13 is completely absent. An important step in this direction is represented 



by the recent paper of ref. [13|, where the problem of mixing has been brought to an 
amenable solution. Here we prove that the calculation can be carried out in a fully 0(a) 
improved way without the need of introducing the Sheikholeslami-Wohlert p8| term in the 



lattice action, or improving the lattice operators [29, 3C] representing the CP-conserving 
effective weak hamiltonian, while still having a positive definite determinant even for non- 
degenerate quarks. 

12 We thank D. Becirevic for drawing our attention on this point. 

13 See, however, ref. p6| for some way-outs if the Wilson theory is fully 0(a) improved a la Symanzik |27], 



4.1 Kaon decay amplitudes 

In the Standard Model the decay of the K meson into pion states is described to leading 
order in the Fermi constant, Gp, by the CP-conserving, AS = 1 effective weak hamiltonian, 



which (in the chiral limit) reads [20] 

"75 



n,AS=l 



V u dV* s — 



C 4 



M w 



M w 



o-(m) 



(4.2) 



where V^KTs is the product of the appropriate matrix elements of the CKM matrix. 14 The 
effective operator TL^ =1 is obtained after having integrated out all degrees of freedom 
above some energy scale, A, with A larger than the charm threshold, but still well below 
the W-boson mass, Jl% Consistently the operators in eq. (|4.2[) are renormalized at the 
scale n with m c <C \x -C Mw, and the Wilson coefficients C±(^/M\y) carry the information 
about the physics between fi and Mw- The expression of the bare operators corresponding 
to 3 ± is 



o- 



{s 1 U){u 1 U)±{s 1 U){u 1 Lu 



1 






u *-* c 


~ 2 





(4.3) 



with 7^ = 7^(1-75). 

For the purpose of making contact with experimental data it is enough to consider the 
decay of the neutral kaon, K°, into either it + tt~ or 7r°7r° states fi~9|, 31, 



the parity invariance of the formal continuum QCD, the relevant amplitudes 



A 



K°- 



A K o_ 



V j=+- 



M w 
M w 



(tt + tt-\O j 
(tt° tt 



VA+AV (H)\K») 



„0| /-/-) j 



Owing to 

(4.4) 
(4.5) 



can be written in terms of matrix elements of renormalized parity odd operators, 
Oy A+AV (fi), with the corresponding bare operators given by 



^VA+AV 



(s7 M u)(u7 / ,7 5 d) ± (s7 M d)(n7 / ,75'u 



u <-> c 



1 








It <- 


-> C 


~ 2 







(4.6) 



4.2 Correlation functions and valence quarks 

Information about the four matrix elements appearing in eqs. (pj) and Q can be ex- 
tracted from suitable spatial Fourier transforms 15 of the correlators (xq > 0, zq > and 

yo <o) 



C 



±,if°7r+7r 



VA+AV 



m J K o (y)) c 



c 



(qcd4) 



Ow0 (x,z,y) = (^o(x)^o(z)0± A+Ay (0)^ (y)) c 



(4-7) 
(4.8) 

14 As usual, the top quark contribution, which is down by a factor 0(VtdV* a JVudVus — 10 -3 ), is neglected. 
15 To get the kaon and the two-pion states at zero (total) three-momentum one has to sum the corre- 
lators ( 4.7)— (4.8) over y and (x + z)/2. In this way the two-pion states can be labeled by the relative 



EE 
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-J 

o 



three-momentum of the two pions. 
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where with the notation (0) CO nn we denote the connected part of the vacuum expectation 
value of the multi-local operator O. As before, the superscript (qcd4) is to remind us that 
correlators are understood to be evaluated in the formal continuum QCD with four quark 
flavours (u, d, s and c). <&k°i ^tt+> an d $ n o are the interpolating operators for kaon 
and pion states, respectively. 

Our approach here is very similar in spirit to that we have proposed in section ||| for 
the evaluation of Br. In order to explain the idea it is convenient to first consider at a 
formal level a gauge model with, as usual, four sea quark flavours (u sea ,, d sea , s sea and c sea ), 
plus several (actually ten) valence quark species. For reasons that will become clear in the 
following, each of the up and charm valence quark species will have to be replicated four 
times. The pattern of valence quarks we wish to consider is then [u, v! , u" , u'" , d, s, c, 
c', c" , d"\. Each quark species is introduced together with the ghost necessary to cancel 
the loop contribution coming from the corresponding fermionic functional integration. We 
will denote the correlators of this model by the superscript (4sl0t>). We now consider the 
(connected) correlators 



C^° v ^_{x,z,y) = ($ n+ (x)%-(z)Q± A+AV (0)& K o(y))< 



(4.9) 
(4.10) 



where the meson interpolating operators $^-0, <& n +, <& n - and & n o are made out of the 
valence quarks u, d and s (but not v! , u" or u'"). This fact will be crucial for the validity 
of the relations ( 4.14| )-( |4.15| ) below. The key point in eqs. d49|)-(^l0|) is the introduction 



of the auxiliary operators Qaja+av^ defined by the formula 



Qt 



VA+AV 



Q ±[0] ,o±[l] 
^■VA+AV + ^VA+AV 



O ^-VA+AV 



±[3] 

VA+AV ■ 



(4.11) 



where the operators Qy jfl^y, k = 0, 1, 2, 3 are all of the form ( [4.6D , but differ from each 
other by the type of u and c valence quark species they contain. Explicitly we have 



^■VA+AV 



>7 M nW)(n[ fc ] 7/ ,75d) ± (57^(^7^75^ 
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. ~ 2 
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l r 



+ 2 



(s^5U lk] )(u [k] ^d) ± (.s~f^d)(u [k] ^u 



1 




~ 2 





(4.12) 



with 



u L ' = u 
cM = c 



u 



[1] 
41] 



u 



uPl = u" . 
A 2 ] = r" . 



= v!" 

J3] = jn _ 



(4.13) 



As argued in the previous section, if the two models (4sl0v) and (qcd4) are regularized in 
the same way (say, using GW fermions), a straightforward application of the Wick theorem 
leads to the equalities 



C 



±,K°n-l 



n _(x,z,y) = C ( ^l +n _(x,z,y) 

n (4,sl0v) , , _ r (qcd4) , , 



(4.14) 
(4.15) 



valid at any finite value of the lattice spacing, provided bare (and thus also renormalized) 
masses of quarks having the same flavour are given the same value in the two theories. 
Eqs. (4.14)-( 4.15| ) in particular imply that the renormalization constants of the operator 



Q-va+av in (4sl0u) can be taken equal to that of the operator Oy A+AV in (qcd4). An anal- 
ogous result is obtained also for the renormalization constants of the meson interpolating 
operators in the two models. Under these renormalization conditions, the equalities ( |4.14 )- 



(4.15) hold for the corresponding renormalized correlators and carry over to the continuum 
limit. 

We conclude this section by offering an intuitive explanation of why we have introduced 
four copies of the u and c valence quarks to get all the mixings canceled. The reason is 
the following. For the purpose of killing the operator mixing coming from "penguin-like 
diagrams" it would have been enough to add to the primary u and c quarks only one extra 



replica of each of them, v! and c', with r u i = —r u and r c / = — r c (see also section 4.5.1 
below). The whole pattern is then duplicated in order to have on each gauge background 
the penguin-like and the non-penguin-like contractions with the correct relative multiplicity 
factor. 

4.3 A convenient lattice discretization of the ci 4 ^ ^ correlators 

±, A 7T7T 

In view of the previous considerations a convenient regularization of the (4sl0f ) model 
is obtained by taking the action fl2.1| ) with valence quarks q u = u, q u > = v! , q u n = u", 
q u m = u'", qd = d, q s = s, q c = c, q c < = c', q c n = c" and q c m = c'" , each having 
an action of the form ([O]) (accompanied by the corresponding ghosts, (f) u , (j) u i, . . . , 4>d, 
4> s , (f> c , 4> c i, ■ ■ ■), plus, as always, the two usual sea quark doublets, ipg = (u sea , d se£l ) and 
iph = (fisea, c sea ) • For reasons which have to do with the need of making as simple as possible 

^VA+AV' 

the valence fields related as follows 



the renormalization pattern of the operators Qy A+A y, we take the Wilson parameters of 



rd = r s , 

r u = -r u > = r u „ = -r u m =r c = -r c , = r c ,i = -r d „ . (4.16) 

For definiteness, we also fix the relation between the two sets of Wilson parameters in 
eq. ( PI) by setting 

r u = r d ■ (4.17) 

The choice r u = —rd is also possible. We comment in appendix [b] on the implications of 
this (isospin breaking) choice. 

The euclidean lattice theory we have just described (see also section ^) will be referred 
to in the following as the (4sl0u)^ m model. 

In order for the correlators of the ^slOu)^ model (with no insertion of sea or ghost 
fields) to be equal to those of (qcd4), we must require the renormalized masses of sea and 
valence fields to be matched as follows 

rhj = rh u = rh u i = m u " = rh u m , rfi^ = rhd , (4-18) 
= m c = m c i = iri c ii = m c m , = rh s . (4-19) 



Owing to eqs. ( |2,6[) -( |2.7| ) and the equality in absolute value of all the valence quark r 
parameters (eq. ( [4.16| )), the formulae (4.18)— ( 4.19| ) fix the non-trivial relations between 



the bare mass parameters of sea and valence quarks and imply the bare masses of valence 
quarks of the same flavour to be equal, i.e. 

m u m = m u ii = m u i = m u , m c /n = m c » = m c > = m c . (4.20) 

As for the interpolating meson fields, one can use the standard pseudo-scalar quark bilinears 
(see, however, appendix |B|) 

*i«> = -^7ss , = -dl5U , = $1+ = uj 5 d , (4.21) 



V2 



(4.22) 



The fundamental reason for considering the (4:slQv)\ m model with its somewhat weird 
valence quark content is that under the conditions ( [4.16| ) and ( |4.20| ) the AS = 1 operators 
Q-VA+AV °f e( l - ( 1^- H| ) m i x neither between themselves nor with any operator of wrong 



chirality or parity. This result will be proved in section 4.5 below by exploiting the symme- 
tries of the (4slQv)t m model. As for the mixing with the famous three-dimensional operator 
sjsd, we will show that, just like in the formal chirally invariant continuum theory, it enters 
with a factor (m^ — m^)(m s — m^), leaving no room for any power divergent behaviour. 
Similarly, despite the breaking of parity at non-zero lattice spacing, no mixing with the 
scalar density operator sd occurs. 

It should also be recalled that in evaluating zero momentum transfer matrix elements 
of Q v a+av 0^- e K ~ > 7F7r )' ^ e °P era t° r sj^d does not play any role, because PCAC allows 
to write s^d = d^s^^^d) / {m s + m^) from which we see that, up to contact terms, s^d 
can always be replaced by a total divergence. 

At this point we can conclude that, after subtraction of the logarithmic divergent 
mixing of Qva+AV W1 ^h (m 2 c — m^)(m s — m^s^d (whenever needed), and multiplication 
by the appropriate renormalization constants, the correlators ( [4.9[ ) and ( f4.10| ), evaluated in 
the (Asl0v)^ m model, admit a well defined continuum limit. The limiting quantities, under 
the conditions fl4.18| )-( [4.19| ), coincide with the (continuum limit of the) correlators (4.7)- 
in the (qcd4) model, thanks to the equalities (fDJ)-(|jJ) . 



4.4 Some observations 

A number of observations are in order at this point. First of all, we stress that, on the 
basis of the general arguments given in section ^, the correlators ( f4.9| )-( ^To| ) can be 0(a) 
improved via W^'s (see eq. ( 2.22[ )). Actually, owing to the symmetry under V x {R — > —R) 



enjoyed by the (AsWv)t m model, the necessary VFj4's of the matrix elements of interest 
can be constructed in terms of data obtained from a single simulation, as explained in 



section |2.3.2| , by averaging — whenever necessary — matrix elements with opposite value 
of the relative three-momentum of the two pions. 

Secondly, for completeness we recall that, in order to actually extract from simulation 
data the amplitudes ■A ft -o_ >7r + 7r - and Axo^o^, it is also necessary to evaluate some extra 



two- and four-point correlation functions, such as 



(<S> K o(x)& KO (x')) , (*')>. <<M*)<M* , )>con n , 

($ 7r+ (x)^ + ( 2 )$ 7r+ (x')< + (2'))coBn , (1> 7r o(x)^o(z)^o(x')^o(z'))conn , (4.23) 

where in the last two four-point correlators one should take xq, zq > and x' Q , z' < 0. Suit- 
able spatial Fourier transforms of these correlators provide the necessary extra information 
on the energies of the kaon, single pion and two-pion states, as well as on the 7T7r-phase 



shift [23, |24[, which will in principle enable us to extract the physical decay amplitudes 
of interest. Finally, when not automatic, the correlators ( 4.23] ) can be 0(a) improved by 
taking appropriate VFA's. 

Once the bare amplitudes of interest have been evaluated with no 0(a) cut-off effects, 
a non-perturbative lattice estimate of the renormalization constants, Zy A+AV (afi), is re- 
quired. The latter can be obtained following the strategy outlined in section 3.3 of ref. Q. 
One could e.g. consider a correlator where QyA+AV * s ms erted at zero three-momentum 
with some other local m.r. operator having appropriate quantum numbers so as not to get 



an identically vanishing two-point function (see also section 2.3.2 ). 

In the same spirit of the observation made at the end of section |3.2| , we note that one 
can eliminate the flavour species u'" and d" from the game (as they are indistinguishable 
from v! and d, respectively) and change the form of the operator in eq. ( [4 . 1 1| ) with the 
substitution 



o± ^o ±[0] +-Q ±[1] --o ±[2] r | o n 

^-VA+AV ^ ^-VA+AV + 2 ^-VA+AV o ^VA+AV ■ N-~ *l 



Again, in this framework the proof of absence of mixing becomes somewhat more compli- 
cated than the one we give below. 

4.5 The renormalization properties of Q-ya+av 

In this section we want to prove that the operators Q-va+av (eq. ( f4.11| )) do not mix with 
any operator of lower or equal dimension, having either wrong chirality or opposite parity, 
precisely as in the formal continuum theory or within chirally invariant lattice regular- 
izations (like those provided by fermions obeying the GW condition). This in practice 
means that no mixing with operators of the same or lower dimension is relevant for the 
computation of the K — ► -kit decay amplitudes. In the arguments of this section a crucial 
role is played by the conditions ( 4. 16| ) and (4.20), which are always assumed as part of 



the specification of the lattice (4sl0v)^ m model we consider. Thanks to eq. ( 4.16 ) and the 
condition |r u | = \rd\, the valence flavour symmetries of the action ( |2.1| ) rule out the mixing 
of Q-va+av an d Q-va+av between themselves (see appendix [C] for details), as in the case of 
untwisted Wilson fermions. We can hence deal "in parallel" with Q va+ av anQl Q-va+av 

4.5.1 Operators of dimension lower than six 

Gauge, H{A) andX/ (with f = s,d, see eq. (fO|)) invariances of the (4^10?;)^, lattice model 
restrict the possible operators of dimension less than six which Qy A+AV can mix with to 



only the following ones 



sd . 



sj 5 a-Fd, 



sa-Fd . 



(4.25) 



where a-F = UnyF^ [U] and F^ U [U] stands for any lattice discretization of the gluon field 
strength. 

We wish to show now that, thanks to the (spurionic) symmetries enjoyed by the 
(4sl0t>)^ m lattice model, the above operators can only come into play multiplied by at 
least three powers of quark masses, according to the pattern given below 



(m c - m u )(m c + m u )(m s - m d ) sj 5 d , 

(m c -m u )(m c + m u )(m s -m d )(m s +m d ) sd, 

(m c - m u )(m c + m u )(m s - m d ) s^ 5 a-Fd, 

(m c - m u )(m c + m u )(m s - m d )(m s + m d ) sa-Fd. 



(4.26) 
(4.27) 
(4.28) 
(4.29) 



As a result, for obvious dimensional reasons none of the above operators can appear mul- 
tiplied by power divergent mixing coefficients. Only the operator ( f4.26| ) will enter with a 



logarithmically divergent coefficient. However, as we have already discussed in section 4.3 



it gives vanishing contributions to the physical decay amplitudes. One might add that the 
mass factor in front of it vanishes anyway in the chiral limit or in the limit of exact SUy (3). 

Since the operators 575^ and s^a-Fd (as well as sd and sa-Fd) have identical transfor- 
mation properties under all the symmetries we will be using in the argument that follows, 
it is sufficient to restrict our attention to the operators (4.2£) and ( [4.27 ). Consistently with 
this observation we notice that the overall quark mass factor in front of the operators ( 4.28j ) 
and ( 4.29| ) is the same as that appearing in front of the operators ( [4.26 ) and ( [4.27| ), re- 
spectively. 

We give a proof of the structure of the mass dependence displayed in eqs. ( 4.26| ) 
and ( f4.27| ) in four steps. 



1. The presence of the factor (m c 



follows from the spurionic invariance of the 



lattice action ( |2.lD of the (AslOv)}^ model under the transformation 



Ex(u,c) x Ex(u',c') x Ex(/,c") x Ex(u'",d") x (m u 



(4.30) 



where Ex(/i, $2) is defined in eq. ( 3.16 ). Indeed, under ( |4.30 ) the operators Qy A+AV 



are odd, while s^§d and sd are both even, as they are left untouched by the transfor- 
mation. 

2. The presence of the factor (m s — m d ), is a consequence of the spurionic invariance of 
the action under the transformation 



1— ' 

o 

O 
O 

O 
-J 

o 



Ex(d, s) x C x (m d <-> m s ) 



(4.31) 



where C is the charge conjugation operation (eq. ( |A.1|) ). Under the transforma- 



tion ( 4.31| ) the operators Q 



± 

VA+AV 



are odd, while sj^d and sd are even. 
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3. At this stage we have proven that both the operators s^d and sd must appear 
multiplied by a factor of the form (m c — m u )(m s — m^). Based on this result we 
now show that a further factor (m c + m u ) must be present owing to the spurionic 
invariance of the action under the transformation 

Ex5(«, v!) x Ex5(V',u"') x (m u — > — m u ) x 
x Ex 5 (c, c') x Ex 5 (c", d") x (m c -> -m c ) , (4.32) 



where Exs(/i, /2) is defined in eq. ( |3.17 ). In fact, under the transformation (|4.32j ) the 



operators Qva+av are even > while the combinations (m c — m u )s^^d and (m c — m u )sd 
are odd. Consequently, the mass factor in front of both operators will have to have 
the form {w? c — m^)(m s — m^). Thus the GIM suppression factor turns out to be 
quadratic just like in the continuum. 

We remark that each of the operators Q-ya+av > k = 0,1, 2, 3, that make up Qy A , AV , 
actually mixes with (m c — m u )(m s — m^s^d and (m c — m u )(m s — m^sd, with linearly 
divergent coefficients. Owing to the spurionic invariance 7£ ' [k] x TZ*?^ (eq. ( 2,16| )), 



such coefficients are, however, odd in r u [ k \ = r c {k\ ■ As a consequence of this prop- 
erty and the peculiar valence flavour structure of the regularization we have intro- 
duced (see eq. ( |4,16D ) , all the mixings with lower dimensional operators cancel against 
each other in the combination ( 4.11| ), which is symmetric under (u,c) <-> (u',d) and 
(u",c") <-» iu'",d"). 

4. On the basis of the further spurionic invariance of the lattice action under the trans- 
formation 

V 5 x (M -» —M) , (4.33) 

where is defined in eq. ( |A.9|) , one concludes that the mixing of Q. VA+AV with sd 
must be suppressed by an extra quark mass factor, besides (mj? — mj) and {m s — md). 
This is so because the operators Qva+av are under the transformation (|4.33j ), 
precisely as (m^ — mj)(m s — m^s^d, while the combination (m^ — mf 1 )(m s — m^sd 
is even. In front of the latter a further mass factor must hence be present, which 
for consistency with the symmetries ( |4.30| ), (4.31) and ( |4.32j ) can only be of the 



form (m s + m^)- We get in this way the expression anticipated in eq. (4.27), which 
we notice is a mere 0(a) cut-off effect. Its contribution actually cancels in Wilson- 
averaged correlators and derived matrix elements. 

We have thus proved that the operators Q-ya+av ^° no ^ m ^ x w ith operators of dimen- 
sion lower than six. 

4.5.2 Operators of dimension six 

The operators QyA+AV ( ec l- (EZC9)) are linear combinations of the operators Qyjf+^y, 
defined in eq. ( [4.12 ). The latter — according to the notation introduced in ( 4.13| ) - 



contain, besides the valence d and s quarks, the flavour species u and c when k = 0, u' and 
d when k = 1, u" and c" when k = 2 and finally u'" and d" when k = 3. 



The somewhat involved proof that the operators Qyjf+^yj k E {0, 1, 2, 3} and, hence 



Q-va+av ( e( l - ( 4.1 1| ) ) mix neither between themselves nor with any other operator of di- 



mension six is given in appendix O. 16 

± \k] 

Incidentally we remark that the overall renormalization constant of each Qy A+A y does 
not depend on the sign of r u [ k ] = r c [ k ] (see eq. Q4.16D ), owing to the general argument given 
below eq. ( |2.21| ). Since all the valence Wilson parameters are taken to be equal in absolute 
value, it follows that the renormalization constant of Qy A+A y is actually independent of 
k, i.e. is the same for all the terms that make up Qy A+A y 

5. The AI = 1/2 rule from K — > n amplitudes 

In the chiral limit K — > tttt amplitudes can be related to K — > tt and K — > vacuum 
amplitudes by the soft pion theorems (SPT's) |33| of Current Algebra [34|. The interest 
of these formulae lies in the fact that they allow to evaluate in the chiral limit K — > tttt 
amplitudes from the knowledge of K — > ir matrix elements which are simpler to compute 
numerically because the final (and the initial) state is a one-particle state and the no- 



go theorem of ref. [22| does not apply. Naturally, as these relations are among amplitudes 
evaluated in the limit of vanishing quark masses, a non-trivial extrapolation to the physical 
mass point will have to be performed. 
Starting again from eq. (f4.2|) 



O that in QCD with four 



it can be proved |1S|, 31, 
flavours it is enough to compute (for a range of quark masses and momenta of the kaon 
and pion states) the matrix elements 



(7T + ( P )\d vv+AA \K + (q)), 



(ir\p)\6 vv+AA \K\ q )) 



(5.1) 



where the bare expression of the parity even operators Oy V+AA reads 

1 r 
2 

+ o (^75«)(«7m75g0 ± {s^^d)(u'~i^u) 



^VV+AA 



(s7 M u)(?x7 M (i) ± (s7„d)(u7„u) 



1 








u <- 


-> c 


~ 2 







+ 



1 








U <r 


-> c 


. ~ 2 







(5.2) 



Following a strategy analogous to the one developed in the previous section, we want 
to prove that the matrix elements ( ^.1[ ) can be extracted (actually with no 0(a) cut-off 
effects) from correlation functions defined in an auxiliary lattice regularized theory where 
the relevant four fermion operators do not mix with any other operator of wrong chirality 
or opposite parity. 

5.1 Correlation functions and valence quarks 

Let us start by noting that the four matrix elements appearing in eq. (fO]) can be obtained 
in the formal continuum QCD theory with four quark flavours (qcd4) from (suitable Fourier 
transforms of) the correlators 



Cgjjfoy) = (^(x)O vv+AA (0)^ K (y)) conQ 

with K and 7r being either K° and tt° or K + and tt + . 

16 As we have already pointed out, the operators QyjfI AV separately do mix with operators of lower 



t 



(5.3) 



dimensions. But in the relevant combination of eq. (4.11) the latter cancel out and arrange themselves so 
as to give rise precisely to the pattern of eqs. (4.26) to (4.29), as was shown in section 4.5.1 



Always at a formal level, we need to consider again the (4sl0v) model that was intro- 
duced in section 4.2. We are interested now in the correlators 



C 



(4sl0u) 
±,Kn 1 



x,y) = ($ n (x)Q 



± 

VV+AA 



mUv)) 



(5.4) 



where the meson interpolating operators §k and 3V are made out exclusively of the valence 
quarks u, d and s, while the operators Qyv+AA are g iven by 



Q 



± 

VV+AA 



^■VV+AA t ^-VV+AA 



o ^VV+AA 



I ±[3] 

^-VV+AA ' 



± [*] 
^■VV+AA 



l r 



u 



[k] 



1 

+ 2 



(s7 M 75« W )(^ W 7 M 75rf) ± (s7M75d)(^ [fcl 7/,75W [fc] ) 





+ 




1 






~ 2 







(5.5) 



(5.6) 



As in the previous section, if the two (4sl0u) and (qcd4) models are regularized in the 
same way (say, using GW fermions), a straightforward application of Wick theorem leads 
to the equalities 

C^ :\x,y) = C^l\x,y), (5.7) 

valid at any finite value of the lattice spacing, provided bare (and thus also renormalized) 
quark masses of the same flavour are given identical numerical values in the two theories. 
Eqs. (|5.7| ) in particular imply that the renormalization constants of the operators Qyv+AA 
in (4sl0v) can be taken equal to those of the operators Oy V+AA in (qcd4). The same 
equality holds for the kaon and pion interpolating fields. With this natural choice of renor- 
malization conditions eqs. ( |5.7| ) are valid for the corresponding renormalized correlators 
and carry over to the continuum limit. 

5.2 A convenient lattice discretization of the C^ 4 ^ ^ correlators 

±, A TV 

For the purpose of getting rid of all mixing problems in the process of computing the 
matrix elements (|5.1|), we will have to regularize the (4sl0f ) model in way which is slightly 
different from that we followed in section |4.3|. Namely, we now have to take 



i— ' 

o 

O 
O 

O 
-J 

o 



-rv = r,,» 



-r s = r d . 
-r u »> = r c ■■ 



-r , = r r ii 



-r c m 



(5.8) 



and fix the relation between the two sets of Wilson parameters in eq. (||~J) by setting 17 

r u = r d . (5.9) 



This lattice regularized model will be denoted in the following as (4sl0u)^. We stress that 
it differs from the lattice model (4sl0v)^ of section |4.3| only by the sign of r s relative to 
r d = r u (and to all the remaining valence Wilson parameters). Naturally, in order to obtain 



17 This choice is made for convenience, i.e. to directly exploit the results of the previous section in the 

exact isospin symmetry in the u and d 



argument we will give in section 5.3 and also to have, for rrid — m 
valence quark sector. However, for the purpose of canceling unwanted mixings only the condition \r u \ = \ra 
is necessary (see appendix H) . 
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from the (4sl0t>)^ model correlators that in the continuum limit are equivalent to those 
of (qcd4), the quark mass renormalization conditions fl4.18|) ~( fi.l9D have to be imposed. As 
in section pP| such conditions entail eq. ( 4.2C] ). 



In practice to extract the matrix elements (|5.1|), besides the three-point correlators 
Q5.4| ), it will be necessary to evaluate in the (4sl0w)^ model also two-point correlators of 
the type 

(^(x)$t( x ')) comi5 {<& K {x)& K {x')) . (5.10) 

Suitable spatial Fourier transforms of these two- and three-point correlation functions con- 
tain all the information which is necessary in order to extract the desired K — » n matrix 
element. The correlators (5.4) and ( 5.10| ) can be 0(a) improved by taking appropriate 



VFA's or by averaging over amplitudes computed with opposite signs of the K and ir 
three- momenta. 

Analogously to the K — > tttt case previously discussed, also here one could eliminate the 
valence quark species u'" and d" and accordingly modify the expression of the operator 
with the substitution 



q± -> Q ± [0] + -0 ± [1] --Q ± [2] mi ) 

^■VV+AA ^ ^-VV+AA ^ 2 VAVAA 2 ^VV+AA • , 



5.3 The renormalization properties of 



VV+AA 
± 



Proving that in the (AslOv)^ model the operator Q VV+AA ( ec l- (HH)) does not mix with any 
other operator is immediate from the results of section [4.5|. It suffices, in fact, to perform 



tin 



in the functional integrals that defines the relevant correlation functions of the (4sl0v) 
model the change of fermionic integration variables induced by the transformation TZ^g (see 
eq. (|2l0|) ). The reason is that under this transformation 1) the action of the (4sl0u)^ 

L 



model turns into that of the ^slOu)^ model except for a sign inversion in front of the 
s quark mass term; 2) the operator Qy V+AA is transformed into the operator Qy A+AV - 
At this point we are exactly in the situation we have discussed in section except that 
we have to carry out the transformation TZ^ S on the list of operators ( 4.26Q to ( f4.29| ) and 



invert the sign of m s in all formulae we previously got. Thus, the pattern of eqs. ( |4.26| ) 
to ( f4.29| ) gets modified as follows 

(m c -m u )(m c + m u )(m s + m d ) sd, (5.12) 

(m c -m u )(m c + m u )(m s + m d )(m s -m d ) sj 5 d , (5.13) 

(m c -m u )(m c + m u )(m s +m d ) sa-Fd, (5.14) 

(m c -m u )(m c + m u )(m s + m d )(m s -m d ) sa-Fj 5 d. (5.15) 



The mass factors appearing in eqs. ( |5,13 ) and ( |5.15| ) in front of the operators 575^ and 



s'jsa-Fd depend quadratically on m s , so they are left unchanged. Furthermore, the whole 
line of arguments we gave before to prove that Qy A+A y could not mix with operators 
of dimension six in the (4-slOv)^ model can be applied to infer the same property for 
the operators Q VV+AA in the (4sl0w)^j model, as dimensionless mixing coefficients are 
independent of the quark mass (in a mass independent renormalization scheme). This 
concludes our proof. 



The net result of this analysis is that the situation is just like in continuum massive 
QCD. Only the particular combination appearing in eq. (5.12) needs a (logarithmically) di- 
vergent coefficient, were this subtraction required. This term, however, will not contribute 
to the K — > ir form factor that is related through SPT's to the K — > irir amplitude, because 
SPT's have to be understood as relations between amplitudes evaluated in the chiral limit. 

Concerning the renormalization constants of the operators Qyv+AA * n ^ ne m °del 
(4sl0t>)^, it is clear that they are equal to those of the operators Qva+av * n ^ ne m °del 
(4sl0t>)^ m (see comments on the evaluation of Zy A+Ay at the end of section O). The 
reason is simply that, as pointed out before, the change of fermionic integration variables 
induced by the transformation 1Z§ S exactly maps the two theories and the two operators 
one into the other. 



6. Concluding remarks 

Our present knowledge of QCD offers a deep understanding of strong interaction physics 
with two noticeable exceptions. One is the lack of a convincing solution of the strong-CP 
problem the second is the difficulty of finding a physically sound explanation of the 
anomalously large value of the K(K — ► irir) (eq. (|4.1| )). 

Clarification of anyone of these two issues may have a significant impact on the struc- 
ture of any forthcoming unified theory comprising the Standard Model, or give hints about 
possible signals of new physics already at low energy, though it is likely that the strong 



CP-problem cannot have a solution within QCD [36|, but requires its embedding in a larger 
theory. 18 

In this perspective the results of the present paper look rather interesting, because, 
when put together with those of refs. ||, |j, set the basis for a workable framework, where 
a first principle computation of the matrix elements of the CP-conserving AS = 1 and 
AS = 2 effective weak hamiltonian at realistic values of the light quark masses appears to 
be feasible with the computer resources that are or will be soon available. A more precise 
estimate of the requested computational effort will be possible as soon as the practical 
issues concerning the magnitude of residual cutoff effects, in particular for light quarks, 
and possible lattice phase transitions in unquenched simulations (see the discussion below) 
will be clarified. 

Independently of these practical issues, the use of maximally twisted Wilson (or rather 
Osterwalder-Seiler) fermions, according to the strategy advocated in this paper, brings 
about conceptual and qualitative progresses which can be summarized as follows: 

1. no zero modes of the Wilson-Dirac operator at non- vanishing quark masses (in par- 
ticular no exceptional configurations); 

2. a positive definite fermionic determinant, even for mass non-degenerate quarks, in 
unquenched studies with an even number of dynamical flavours (in this paper we 
focused on the phenomenologically relevant case of four flavours); 



3 See, however, the papers in ref. [^] for a possible solution within QCD itself. 



3. no contribution from wrong chirality or opposite parity mixing for the lattice opera- 
tors entering the computation of the weak matrix elements of interest; 



4. automatic (or easily obtainable) 0(a) improvement of all relevant physical quantities. 

Although the method we propose is rather general, in this paper for definiteness it 
has been applied only to the (important) instance of the evaluation of the Bk parameter 
entering the K°-K° mixing and the K — ► tttt, or K — > n, amplitude. In the lattice 
computational framework we propose 

• the Bk parameter can be hopefully computed including sea quark effects, with a 
sufficiently small error, so as to significantly reduce one of the largest theoretical 



uncertainty in the phenomenology of the unitarity triangle [15 



• the yet elusive AI = 1/2 signal does not appear to be anymore buried under- 
neath a numerically overwhelming, divergent subtraction (as was the case for Wilson 
fermions) . 

These results are obtained at a surprisingly low price. Sea quarks need to be introduced in 
pairs and a somewhat exotic pattern of valence quarks has to be used in the construction 
of the effective weak hamiltonian lattice operator. 

From several viewpoints the situation is thus very much like with GW fermions, where 
chiral symmetry is exact and cut-off effects only show up at 0(a 2 ). The non- negligible gain 
that we get with the strategy we are advocating here is that the computational burden ap- 
pears to be significantly reduced with respect to simulations involving GW fermions. This 
feature makes our method potentially suited for the study of systems with a physical linear 
size of a few fermi, while including sea quark effects. In comparison to approaches based on 
untwisted Wilson fermions (which are, however, plagued with complicated operator mixing 
and the problem of "exceptional configurations" ) , the method we propose possibly requires 
a little more extra computing time, owing to the need of computing certain valence quark 
propagators for opposite values of the Wilson parameter. This overhead is more than com- 
pensated at small quark masses by the expected improvement in the performance of linear 
solvers, due to the protection against spurious quark zero modes that is guaranteed by 
twisting. 

6.1 Critical mass and O(a) improvement 

It should be noted that our work is not based on any special definition of the critical 
mass for sea and valence quarks (see eqs. (|2.4[)-(|2l)|)). The method for 0(a) improvement 



discussed in section 2.3, which is a straightforward generalization of the approach presented 
in ref. H, relies only on the property that f C r( r i] r 2> r 3) is an odd function of r%, implying 
that the sea and valence critical masses, M cr (r^h) and M cr (rj) are in turn odd in and 



Tf, respectively. The argument showing that f cr is odd in r\ we give here (in section 2A 
and appendix |A]) is analogous to that presented in ref. [g] in the case of generic (chirally 
twisted) Wilson fermions. 



The question of the r-parity properties of M cr (r) was taken up in refs. [38, 39]. Clearly 
this issue could have an impact on the question of 0(a) improvement as we have discussed 
it in the present paper and in ref. ||]. Actually, it turns out that M cr (r) can and should 
always be taken odd in r. 

The reason is simply that any counterterm to be included in the process of renor- 
malizing the theory should not break any of the symmetries the classical action enjoys 
prior to renormalization. Indeed, the manifold of the possible solutions for M cr , defined 
by whichever condition one wishes to take, has the property of being reflected around zero 



under the transformation r — > —r (see e.g. ref. [4C] for a discussion of this point). Thus it 
would be unwise and not very useful to violate the (spurionic) symmetries of the theory 
(in particular Px(i!-t — R)) by having a critical mass that, for instance, to some 0(a p ) 
with p > 1, is not odd in r. 

Apart from this observation, it is important to remark that the whole issue of 0(a) 
improvement, as it is discussed in ref. ||] and in this paper, assumes that as a — > the 
counting of powers of a is done by taking the renormalized quark masses fixed. We did 
not address the interesting question of what happens when the quark mass is taken to be a 
quantity of 0(a) or 0(a 2 ), as often done in papers dealing with lattice chiral perturbation 
theory |§ B S © @- 

6.2 Tests and numerical issues 

The numerical effectiveness of the whole approach based on maximally twisted Wilson 
quarks has to be carefully checked, especially for small quark masses and when taking into 
account the effects of sea quarks. 



There are already encouraging results, coming from the quenched simulations of ref. [12] 
and the analytic work of ref. pl| , indicating that the scaling behaviour of the theory is 
indeed smoother and flatter than with standard or clover-improved Wilson fermions. An 
extension of the scaling tests of ref. |l2| to several lower values of the pseudoscalar meson 



mass, down to less than 300 Me V, is currently in progress [44|. 

Furthermore, recently, an exploratory study of the effects of unquenching in simulations 
of tm-LQCD with two mass degenerate quarks on rather coarse lattices of size 8 3 x 16, 12 3 x 
24 and 16 3 x 32 at j3 = 5.2 (using the standard Wilson plaquette action) was presented Eq|. 



In this work the PCAC quark mass, m^ CAC (see eq. (12) of ref. EH]), the mass of the pion 



and the plaquette expectation value have been measured at ami = 0.01, 19 as functions of 
the untwisted mass, mo. 20 



Results show features of the kind suggested by refs. J46], 42, [43| in the framework of 
chiral perturbation theory, when the order of magnitude inequality a 2 AQ CD <C amp is not 
satisfied Q. Metastabilities are detected, signalled by the behaviour of the plaquette expec- 
tation value which is seen to undergo hysteresis cycles as tuq is driven through M CT . In this 
situation gauge configurations cannot be considered as properly thermalized. Metastable 
states associated with values of m^ CAC of different sign have been identified and argued to 

19 In the notation of ref. m< = fi. 

20 Introducing mo in the action ( |2.l[ ) amounts to replacing there W C r with W„ — M cr + mo. 



stem from the effective lack of criticality of the subtracted Wilson term and the associated 
existence of non-zero values for the parity and flavour breaking quark condensates. As a 
consequence, one does not yet observe the expected vanishing of m^ CAC at a critical value 
of mo (providing a natural definition of M cr ). It is, however, reassuring to see that the 
squared pion mass appears to be roughly proportional to J (ZAm^ CAC ) 2 + mj in both the 
observed metastable branches. 



On the basis of the present understanding [43| of the metastability phenomena dis- 
cussed above, we expect them to get weaker and weaker and eventually disappear as the 
lattice spacing is decreased, making thus possible the evaluation of physical quantities with 



controlled 0(a ) cut-off effects. Preliminary data [47] coming from simulations in progress 



at /3 = 5.3 seem to be consistent with such an expectation. 
6.3 Outlook 

The idea of using a hybrid approach where valence and sea quarks are regularized in 



different fashion was already put forward in the literature (see, for instance, ref. |4S] and 
references therein). Along this line a winning strategy, where all the nice properties we have 
listed above remain valid, could be to use GW fermions as valence quarks and maximally 
twisted quark pairs in loops. With this choice one could ease the problem related to the 
need of working at relatively large values of (3 (and correspondingly large lattices) if the 
valence quark mass is small, and obtain a computational framework, where CPU times 
for inverting the overlap Dirac operator and stochastically taking into account the fermion 
determinant are comparable. 

We end by recalling that the tm-LQCD action can also be employed Q in the frame- 
work of the Schrodinger functional formalism |50[| . Results about the 0(a) improvement 
remain valid for quantities that are independent of the boundary conditions, like on-shell 
matrix elements or masses. However, if one wishes to improve the whole Schrodinger func- 
tional, which is a necessary step in order to get cancellation of 0(a) cut-off effects in the 
evaluation of, for instance, the running gauge coupling, quark masses and renormalization 
constants, the way of extending the WA method is not completely obvious and will be the 



object of a forthcoming publication [51|. 
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A. Renormalizability of the lattice model ( |2.1|) 



We want to show in this appendix that the action of the lattice model (|2,l| ) contains 
all the terms of dimension not larger than four, allowed by the action symmetries, that 
are required for correlators with no ghost fields be renormalizable. In the course of this 
discussion we also prove 1) the results contained in eqs. ( |2.6[ )-(2,7), 2) the property that 



a unique function, / cr (r; r/>, r^), determines the critical mass of OS valence and twisted sea 
quarks (see eqs. (|2.4j) - (|2.5[ )) and 3) the relation Z m {rj) = Zp[rf) for 77 = r£ or rj = r^. 
We start by noticing that completely standard arguments, based on dimensionality 



and invariance of the action (2.1) under lattice gauge transformations, translations, H(4) 



rotations and charge conjugation, 

" u,(x) - u;(x) 



C : < 



t/) p (x) i~fol24>p(x) , i> p {x) -> -^J(x)i7o72 , p = £,h 
q f (x) -> ^70729/ (x) , qf(x) -> -qj {x)ij -f 2 , all /'s 
. 4>f{x) i7o72(0/) T (2;) = «7o72^}(x) , all /'s 



(A.l) 



imply that the only relevant action density terms are, besides gauge invariant gluonic 
operators (FF and FF), either fermionic bilinears of the type i/L . . .ipp and qt . . .qf or 
ghost bilinears of the form <jy. . . .<f)f. 

A.l Sea quark sector 

The renormalization of correlation functions with only gluons and ijj£,h fields proceeds 
as discussed in ref. In fact, for this class of correlation functions the contributions com- 
ing from functionally integrating valence quark and ghost fields completely cancel against 
each other, while the symmetries discussed in ref. are straightforwardly extended from 
one to two pairs of maximally twisted quarks. 21 

Furthermore mixed sea quark action terms of the form ipi . . . tph or iph . . . ipg can not 
arise, as it trivially follows from the invariance of the action under the independent U(l) 
vector transformations 

X v : il> v {x) -» e ie *i>p(x) , $ p (x) -► ^(x^"^ , p = l,h. (A.2) 



The important equation ( |2.q ) can be derived from the analysis of the chiral WTI's 
of the lattice model fl2.1| ) involving correlators with the insertion of ij)i h and ^ ^ fields 
only. Since valence quarks are totally inert, the structure of these WTI's is identical to 
that reported in ref. Q for maximally twisted lattice QCD with one pair of mass non- 
degenerate quarks. For the reader convenience we report here the WTI's from which the 
renormalization properties of the quark masses most straightforwardly follow. If 0(y) is a 
renormalized (multi-local) operator, one gets, among others, the two WTI's (x ^ y) 



(d;Af(x) - 2m p P^{x))d{y)) |^ = 0(a) , (A.3) 
d*A 3p (x) - 2m p P 3p (x) + e p P 0p {x))6{y) 



(R,M) 

= 0(a) , (A.4) 



(R,M) 



21 In particular the symmetry operations that do not leave inert the gauge field must be performed on 
both quark pairs simultaneously. Important examples of these symmetries are 

Vl/2 x (mj -* -mi) x (m h -> -m h ) , V% x (ej — > -ej) x (e h — > -e h ) , 

with V\ 12 and V% representing the obvious extensions of the operations in eqs. (62)-(63) of ref. ^| to the 
theory with two quark flavour pairs. 



where p runs over all possible sea quark pairs present in the theory, p = £,h. In eqs. ( |A.3| )- 
( |A,4| ) we have used the definitions 

Af? = Z v Vyy^y V>p , Aff = Z v Vv7m75tt , (A.5) 
P 3 p = Z P Vv^ 7 5^ P , P 0p = Z s ^75^ , (A.6) 

and d* is the backward derivative. The above formulae imply the relations 

rhp = rhp + ip = Zp^m v + Z<7 1 e p , (A. 7) 

rh p = rhp — i p = Zp^m p — Z~^e p , (A. 8) 

for p = I and p = h. 

The renormalization constants Zy, Zp and Zg are even functions of the Wilson pa- 
rameter of the quark pair we focus on and the r parameters of all the other sea quark pairs 
appearing in virtual loops, Zp[r p ) = Zp[r p ;ri,rh), I = V, P, S, ... The renormalization 
constants Zp[r p ) coincide with the corresponding ones of the standard Wilson theory with 
four flavours, if a mass independent renormalization scheme is employed. 

A. 2 Valence quark sector 

Let us now come to the proof that the valence quark sector of the action ( |2.1| ) has the most 
general renormalizable form compatible with the symmetries of the model. 

Due to valence flavour conservation, each valence flavour / renormalizes independently 
from all the others. Moreover one verifies the validity of the following statements. 

1. The d = 4 operator g/757 • Vqf is forbidden by C; 

2. The d = 3 operator iqfjsqf can appear with a coefficient (to be identified with the 
critical mass), proportional to 1/a and odd in rj. Such a term is allowed by the 
spurionic invariance V5 x (M — > — M), where 



^5 : < 



U {x) -» U (xp), U k (x) -» Ut(xp - ak) , k = 1,2,3 

t/j p (x) -> j57o^p(x P ) , ip p (x) -> -Vv(xp)7o75 , p = £,h 

q f (x) -> 75709/ (xp) , q~f(x) -> -g/(a;p)7o75 , all f's 

, 4>f(x) 75700/ (xp) , all f's , 



(A.9) 



as such symmetry only excludes coefficients odd inrnj. The spurionic invariance IZ^ 
(see eq. ( 2 . 1 6[ ) ) fixes the rj-parity of the coefficient to be negative. We recall that no 



odd powers of mass can arise from the sea quark determinant, as the latter is even in 
the sea quark mass parameters (as well as in the sea Wilson parameters). This proves 
the statement, made in the text, that the critical mass M cr (rj) = a" 1 / cr ( r f\ r ii r h) is 
an odd function of Tf. 

3. The d = 3 operator cjfqf is allowed, owing to the spurionic invariance X (M — > 
—M), only if multiplied by an odd power of mj (again no odd power of mp^ or 
can arise from the sea quark determinant). Hence a term like a~ 1 qjqj is ruled out 



and only the term mfqtqf needs be included in the action. This implies the multi- 
plicative renormalizability of the valence quark mass parameter, m/, see eq. ([2.7|). 22 
As remarked before, the renormalization constant of m/, called Z m {rf) in the text, 
must be even in rj, for consistency with the spurionic invariance Tx^j- 

4. The reality properties of the coefficients of all the quark bilinears are fixed by site/link 
reflection invariance, O s //. The action of Q s /i on gauge and sea quark fields has been 
specified in eqs. (58)-(59) of ref. [f§. Valence quarks, qf, and ghosts, 4>f, trans- 
form under Q s ji precisely as each member of the sea quark pair field, ip p , while qj 
transforms as the corresponding antiquark field. 

Since we are not interested in correlators with ghost field insertions, we need not discuss 
in detail the renormalization of the ghost sector of the model. We only wish to stress 
again that, as we already remarked in the text, the parameters, rj and rrif, in the ghost 
action fl2.3| ) are completely fixed by the requirement that integration over the ghost field 
should exactly cancel the valence fermion determinant. 23 

It is also a trivial observation that the invariance of the model ( |2.1| ) under independent 
U(l) transformations acting either on single valence flavours (If, see eq. ( |2.9D ) or individual 
pairs of sea quarks (I p , see eq. ( |A.2| )) guarantees that no fermion bilinears with mixed 
valence-sea quark content can arise. 

As for the operator iFF, its presence can be induced neither by integrating over the 
sea sector Q, because of the invariance under parity (V, see eq. ( 2.26| )) combined with 



4>l,h{x) -> iT 3 ip i)h (x) , ipi,h(x) -> -tp£,h(x)ir 3 , 

nor by integrating over the valence sector because the resulting fermionic determinant is 
exactly canceled by the ghost integration. 

We note incidentally that the transformation V5 (see eq. (A.9)) can be also cast in the 
form 

V h = V x K u x TZ 5h x JJ n sf , (A.10) 

/ 



where the transformations TZ^f and TZ^ P are defined in eqs. ( 2.1C| )-( |2.11| ). 



A. 3 Sea and valence quark mass renormalization 

Here we want to establish the fact that within the model (2.1) a unique dimensionless 



function, / cr , determines the critical masses in eqs. (p.4|) — ( p~5| ) and prove the relation ( |2.§| ) 
between the quark mass renormalization constants of sea and valence quarks. 



22 The analysis of chiral WTI's with insertions of valence quark operators leads to the same result. 

23 Incidentally if one is willing to consider graded field transformations between valence fermions and 
ghosts [jnj, then it turns out that, once the valence fermion action is fixed by the usual renormalization 
conditions, the form of the ghost action is determined by requiring invariance under e.g. the transformation 
(with super-determinant equal to 1) <j>;(x) — > qf(x), qf(x) — * 4>f(x), 4>\{x) —* sign(m/)g/(x), q~f(x) — > 
sign(m/)(/>J.(a;), which for each flavour exchanges quark and ghost fields. 



For this purpose, let us consider correlators containing either only the fields ipi and 
or only a pair of valence flavour fields, that we collect in the doublet Q = (01,(72), an d 
choose 

r i = r i = ~ r 2 = r j m £ = m% = mi = m , eg = . (A. 11) 

In this setting the relation between quark mass renormalization constants we wish to prove 
reads 

Z m (r;r,r h ) = Zp X {r\r,r h ) . (A. 12) 

It is convenient to consider two pairs of WTI's, one involving correlators made out of 
exclusively the fields ipg and ip£ and the second one involving only the valence fields Q and 
Q. The WTI's of interest are (x ^ y) 



d;(Al*(x)P 2e (y)) 



(R,M) 



(R,M) 

where we have used the definitions 

v™(x) = Mx)^Mx), 



0(a), 

2m(P 2e (x)P 2e (y)) 



(R,M) 



+ 0(a), 



p 2e (x) = Mxhs^M*) , 



T2 



2A l e ( x ) = i'l { x )ltil5 Ufi (x)ij)£ (x + a/2) + ij)t(x + 0^)7^75— C/^(a;)^(x) 

_ 7~3 — T3 x 

^t^ — U^x)^^ + a/2) - ifj e (x + ap,)—Uj l (x)il)e(x) 



(A.13) 
(A.14) 

(A.15) 
(A.16) 

(A.17) 



and the corresponding WTI's with the ipg and ipg fields replaced by the valence fields Q 
and Q, i.e. 



d;(V 2Q (x)P 2( 3(y)) 
dl(A 2 ^(x)P 2 ^(y)) 
where analogously we have set 



(R,M) 
(R,M) 



0(a), 

2m(P 2Q (x)P 2Q (y)) 



(R,M) 



+ 0(o) 



T2 



P 2 ®(x) = Q( x ) lb l±Q( x ) , 



T 2 . 



2 4«(x) = Q(x)7/,75y^(x)g(x + a/i) +Q(x + aA)7 M 75y^(x)g(x) + 



T-2- 



+r 



Q{x)^{x)Q{x + a/2) - Q(z + ap)^Vl{x)Q{x) 



(A.18) 
(A.19) 

(A.20) 
(A.21) 

(A.22) 



Imposing the validity of the WTI's ( A.13| ) and ( A.18| ) is a way of arriving at a definition 
of the critical mass for the fight quark pair and the valence quark doublet Q = (01,02), 
respectively. The key observation is that, with the choice ( |A.ll ) of bare parameters, the 
l.h.s. of eq. ( A.13| ) is exactly equal to the l.h.s. of eq. ( A.18| ) . To prove this it is enough 



to perform a (suitable) vector- T2 rotation of the fight quark field pair, so as to rewrite its 
action in the basis where the Wilson term is multiplied by —175X3 as is the case of the Q- 
doublet. Such a rotation does not alter the expression of the operators ( A,15| ) and ( |A,16j ), 
but brings the flavour structure of the axial current ( |A.17| ) to that of eq. ( |A.22| ), It thus 
follows that the numerical value of M cr = a _1 / cr (r; r, r^) inferred by enforcing the vanishing 
of the r.h.s. of eq. ( A.13| ) is the same that is obtained by enforcing the same condition on 
eq. ( [A. 18 ), and viceversa. This proves the relations (2.4)-( |2~5| ). 

The WTI's (A. 14) and ( A.19| ) contain all the relevant information about the quark mass 
of the fight quark pair and the valence quark doublet (91, 92), respectively. Moreover they 
involve the exactly conserved 1-point split axial currents ( ]A.17 ) and (A. 22). Precisely for 
this reason, the quark mass renormalization factor can be taken equal to the inverse of the 
renormalization constant of the relevant pseudo-scalar densities, which are P 2i (eq. ( A.16| )) 
and P 2Q (eq. ( |A.2l|) ) for the fight quark pair and the valence quark doublet, respectively. 
In the text we have denoted the mass renormalization constant of the fight quark pair by 
) = Zp^ir^r^Th) and that of the valence fields Q = (91,92) by Z m (r) = Z m (r;r, 77J 



Zl 



(recall that Z m is even in r and we have taken r± = —V2 = r, see eq. ( A.ll| )). At this point 
the relation ( |A.12 ) is a consequence of the fact that, by the same argument we gave before 
for the vector WTI's ( gjg ) and flATs)) , the correlator (P 2i (x)P 2i (y)) \ {RM) in eq. ( |Al4l ) 
is equal to the correlator {P 2 ®(x)P 2 Q(y))\( RM ) in eq. ( |A.19| ). 

B. Isospin and parity breaking effects in the light quark sector 

Although to guarantee the desired renormalization properties of Qy^^y only the condi- 
|r^| is actually necessary, the choice (4.17) is quite convenient as it preserves 



tion \r„ 



(for rrid = rn u ) the isospin symmetry of the action of the u and d pair of valence quarks. 
As a consequence of this symmetry, the two-pion states \-k + tt^ > and 1 7r°7r° > cannot mix 
with the neutral pion state, \ir° >, as it would be the case had we chosen the condition 



r u = —r<i instead of (4.17). This remark is very important in practice, because (recall that 
parity is not a symmetry of tm-LQCD) it allows to exclude the presence of the neutral 
pion state among the low mass states contributing to the spectral representation of the 
correlators C^+^_ (x, z, y) and C^^o z, V) (see eqs. (p|)-(|4.iq)). 

It should be noted, in fact, that, if present (as in the case isospin is broken by the 
choice r u = —r^ or by taking ^ m u ), the contribution of the neutral pion may numeri- 
cally compete with the relevant isospin and parity preserving contributions, since its mass 
is generically lower than the energy of the two pion states we are interested in. Consis- 
tently with the general properties of tm-LCDQ it can, however, be shown that isospin and 
parity breaking contributions to the spectral representation are mere 0(a 2 ) cut-off effects. 
The latter can be disentangled by standard techniques and removed from the measured 
correlators before extracting the relevant 0(1) contributions. 

There is an extra problem that should be mentioned if isospin (and parity) is broken. 
Namely the fact that the standard interpolating field for the neutral pion, eq. fl4,22j ), mixes 
with the identity operator. To avoid this problem one may replace the latter with the 
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fourth component of the neutral axial current 



1 

71 



7/707514 - d7o75^ 



C. Absence of the mixing of Q VA+AV with operators of dimension six 

In this appendix we prove that the operators Qy^+^y cannot mix with any other operator 
of dimension six. Looking at eqs. ( f4.16| )^(4.17), we see that just two operators need to be 
considered, i.e. 



1- QyA+AV > wnere valence quark species have r u 



2. QyA+^yi where valence quark species have r u > 



r r > 



-rd 



"TV 



Indeed the case of QyA+AV * s identical to the case of Qyjf+^y, an< ^ ^ ne case °^ 

to that of Qyjji^y, because the quarks that enter have an equal pattern of signs of the 

Wilson parameters. 

Since we are discussing the mixing of operators of the same dimension, we can set to 
zero the masses of the quarks and consider the lattice model {AslQv)\ m in the chiral limit. 
Let us discuss separately the case of Q^i^y and Qyjjl^y. 



C.l The operator QyjJ^jy 

A basis for the operators of dimension six which Qyji+^y can mix with is given by 

QP 0] , Qy [0 \ 

Y = VA + AV, VA - AV, SP + PS, SP - PS, TT, 
VV + AA, VV-AA, SS + PP, SS-PP, TT . 



(C.l) 



The notation VA + AV, VA — AV, . . . , TT refers to the Dirac structure of the four-fermion 
operators and is self-explanatory. For instance, with a^ u = i[y^, , j u ]/2, we have 



Q±[0] 

^■VV+AA 



1 



+ - 



(s7 /1 u)(u7 At d) ± (s7 At d)(ti7 At n) 



1 



\u c + 



(s7 jU 7 5 'u)(ti7 At 75<i) ± (s7^75^)(^75^) 



u 



Q 



±[0] 
TT 



— \u <-> c . 
2 L J 



(C.2) 
(C.3) 



The key observation is that, owing to the equality r u = r c = = r s of the Wilson 
parameters of the quarks appearing in Qy'°', the part of the lattice action relative to 
these quarks (the standard u, d, s and c quarks) admits an exact SU(4) flavour symmetry, 



The arguments that follow would go through also if we were to replace eq. ( ^.lTf ) with r u = — r<j, while 
maintaining eq. (4.16). In fact, this new choice for the valence Wilson parameters is related to the one 
we made in the text by the transformation Ex(u, u) x Ex(tt",u'") x Ex(c, c') x Ex(c",c"') which leaves 



invariant the form of the operator Q 



VA+AV 



See, however, appendix 



which prevents the operators of the type (+) to mix with any of the (— ) counterparts, 
as they belong to different irreducible representations of the SU(4) flavour group. The 



operators in (CA_) can be thus viewed to yield two separate bases, one for the (+) sector 
and another one for the (— ) sector, containing the candidates for the mixing with Qyjy^y 

and Q V a+av> respectively. 

Notice that conservation of each separate valence quark species forbids the mixing of 
the operators Qyjj'i^y with the operators Qy^ k \ k = 1,2, 3, as they have different flavour 
content. 25 

For brevity, we will deal "in parallel" with Qyj^+^y and Qyjy^y One checks that 
the invariance of the action (&slQv)\ m under V5 (see eq. ( |A.9| )) rules out the mixing of 



QyA+AV w ith the operators of the form Q^} u> , having the Dirac structure W = VV + 
AA, VV — AA, SS + PP, SS — PP, TT. In fact, these operators are even under V5, while 

Qva+av is odd - 

To complete the proof we notice that under the symmetry transformation Ex(d, s)xC, 
the operator Qyjj^y is odd, while the operators Q^p+ PS and Q^~°^ are even. As for 

the operators <2yjy_^y and Q^p}_ PS they have more complicated transformation proper- 
ties, namely Qyjjl^y — ► Qyjy^y and Q^p^_ PS — ► — Qjp°l P5 . The operators that under 
Ex(<i, s) x C have the same transformation properties as Qyjf+^y are thus the combina- 
tions QvaLav — Qyjy^y and <2^p°l P5 + Q SP °}_p S , but SU(4) flavour symmetry prevents 
Qyjl+^y from mixing with them. This concludes the proof for Qyjy^y 

C.2 The operator Qy [ A 1] +AV 

We now consider the case of the operator Qyj^+^y, which differs from the previous one 
only because the Wilson parameters of the relevant valence quarks obey the relation r u > = 
-r d = -r s = r d . 

The mixing and renormalization properties of the operator Qyju^y can be brought 
back to that of the operator Qyj^i^y, discussed above, by simply observing that the change 
of functional integration variables induced by the transformation TZ^ U ' x TZ^ does not touch 
the expression of Qyji+^y and leaves invariant the massless action of the valence quarks, 
apart from switching the signs of the Wilson terms of the quark species v! and c'. We get in 
this way precisely the situation we had before. Since a change of variable in the functional 
integral leaves unchanged all the correlators, from the absence of mixing we proved in the 
case of Qyj^^y we can infer the same result for Qyy^y- This ends the proof. 
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